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We have calculated the absorption characteristics of different hybrid systems consisting of Ag, Ag 2 or Ag 3 
atomic clusters and poly(methacrylic acid) (PMAA) using the time-dependent density-functional theory. The 
polymer is found to have an extensive structural-dependency on the spectral patterns of the hybrid systems 
relative to the bare clusters. The absorption spectrum can be 'tuned' to the visible range for hybrid systems 
with an odd number of electrons per silver cluster, whereas for hybrid systems comprising an even number of 
electrons, the leading absorption edge can be shifted up to ~ 4.5 eV. The results give theoretical support to the 
experimental observations on the absorption in the visible range in metal cluster-polymer hybrid structures. 
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I. INTRODUCTION 

When metal samples are reduced from bulk to nanome- 
ter dimensions, their optical properties become governed 
by surface plasmons. When the size is further decreased 
to few-atom clusters, the properties are again altered as 
discrete energy levels determine the electronic and opti- 
cal behavior of the cluster 1,2 . Applications in this size 
regime have become extremely important during recent 
years as advanced experimental production on this scale 
has become more controlled and theoretical methods are 
now more easily able to access it 3 ' 4 . 

Silver clusters are especially interesting but challeng- 
ing materials. Their high biocompatibility suggests great 
potential as functional components in biological nanoma- 
terials with possible applications in biosensing and op- 
tical data processing '~ 10 . Silver clusters readily oxidize, 
which can be a problem for some applications, and meth- 
ods of increasing their stability against oxidation are ac- 
tively being sought 11 . From the theoretical perspective, 
the challenge of these systems lies in the hlled d-electron 
levels of silver, which makes high requirements for model- 
ing their electronic properties, and for determining their 
potential for application. 

During the last two decades, both extensive experimen- 
tal and computational studies have been done on small 
silver clusters (Ag„) with less than a dozen atoms. Ex- 
perimentally these clusters have been studied in stabi- 
lizing rare-gas matrices, as free silver clusters in the gas 
phase dissociate by fragmentation 12-1 ^. Small silver clus- 
ters exhibit interesting luminescence and especially fluo- 
rescence properties that can be tuned by their size and 
chemical environment 16-1 8 . Interesting, but not very well 
understood, experimental observations on the tunability 



of their optical properties as a function of the chemical 
environment have been reported. Di'ez et al, for example, 
have observed absorption in the visible range at ^500 nm 
(2.5 eV) from small silver clusters (n = 1 — 4) stabilized 
by coils of poly(methacrylic acid) (PMAA) in the liquid 
phase . The photoabsorption spectrum of this system 
was found to rcdshift slightly when the polarity of the sol- 
vent decreased, and other interesting properties, such as 
strong fluorescence and electrochemiluminescence, were 
also observed. A similar synthesis was also carried out 
by Shang et al. 20 who assumed that the measured ab- 
sorption bands at 430 nm (2.9 eV) and 500 nm (2.5 eV) 
are due to Ag+ clusters with n = 2 — 8. Their study also 
reported strong fluorescence. 

On the computational side, methods used to model 
excited-state properties of nanosized structures include 
the coupled cluster method 21 , ab initio configuration 
interaction (CI) approach--, time-dependent density- 
functional theory (TDDFT) 23 , Hedin's GW approach 
combined with the solution of the Bcthc-Salpeter equa- 
tion 24 , and the quantum Monte Carlo method 2 '. The 
quantum chemistry community, in particular, Bonacic- 
Koutccky et al. pioneered work on the absorption proper- 
ties of silver clusters using ab initio CI, publishing many 
bench-mark results 22 ' 26 ' 27 . On the DFT side, Yabana 
and Bcrtsch were the first to use a real-time implementa- 
tion of TDDFT to calculate the optical response of small 
silver clusters Ag 1; Ag 2 , Ag 3 , Ag 8 and Agg in the gas 
phase 2 *. A more systematic real-time TDDFT study of 
the optical properties of Ag n with n = 1 — 8 followed 21 . 
Both of these studies showed good agreement with exper- 
imental results with a rather rigid difference of about 10% 
in the excitation energies. More recent TDDFT studies 
have been performed on larger Ag n clusters with n > 4 
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A hybrid structure consisting of a silver cluster and an 
organic compound is of great interest as this is already a 
reliable model for a silver cluster that has been function- 
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alizcd by its environment. Examples of studies on silver- 
cluster-organic-compound hybrid systems include silver- 
carboxylates for which DFT methods have been used to 
determine the shapes of the resulting silver nanoparti- 
cles 33,34 . The interaction with the organic compound 
can cause strong modifications to the optical properties 
of the cluster. For example, Mitric et al. studied the ab- 
sorption spectra of silver cluster-tryptophan hybrid sys- 
tems using the TDDFT 3 ". They observed a strong size- 
and structural-dependence of the photoabsorption and 
photofragmentation in these systems. Other examples of 
recent TDDFT work on silver-cluster-organic-compound 
hybrid structures include a study of the enhancement 
mechanism in surface enhanced Raman scattering 36-38 
and the determination of the optical properties of silver- 
clustcr-biomolccule- hybrids . 

In this paper, we study the effect of organic compounds 
on the photoabsorption spectra of silver cluster systems. 
For our model systems, we choose clusters consisting of 
one to three silver atoms attached to fragments of the 
PMAA polymer. This choice is motivated by recent ex- 
periments in Refs. 19 ' 20 ' 40 where small silver clusters have 
been formed in the presence of PMAA in different so- 
lutions. In our calculations we focus on the cluster- 
polymer system only and omit the solvent for simplic- 
ity. This approach has also been used by Mitric et al. 
who calculated the excitation properties of interacting 
nanoparticlc-biomolecule subunits . 

The paper is organised as follows. The selected ge- 
ometries for our systems are investigated and presented 
in Sec. II. The DFT and TDDFT methods are covered 
in Sec. III. Sec. IV contains the results and the analysis 
of the ground-state calculations and Sec. V contains the 
results of the optical absorption (excited-state) calcula- 
tions. Finally, the conclusions are given in Sec. VI. 



II. STRUCTURES 

In order to emphasize their essential features, a 
schematic representation of the structures we have stud- 
ied is shown in Fig. 1. Full geometry optimization on 
these systems have also been performed, and the de- 
tails of these calculations can be found in Sec. III. The 
acronyms used for these structures consist of a descrip- 
tive letter (P = polymer, M = monomer Agi (atom), 
D = dimer Ag2, T = trimer Ag3) followed by a con- 
secutive number. Our choice of systems for the hybrid 
structures is motivated by experimental evidence, which 
shows that carboxylate groups are capable of bonding 
with Ag ions in polyacids (see, e.g., Ref. 4 " and references 
therein). Hence, in this work, our primary interest has 
been to study small (i.e., single atom and dimer) hybrids. 
For the monomer (M1-M8) and dimer (D1-D7) struc- 
tures, a large variety of different geometries have been 
systematically studied, whereas only a few trimer struc- 
tures (T1-T2) have been used to test whether the trends 
obtained for monomer atom and dimer structures can be 



extended for systems containing larger Ag n clusters. 

The pure polymer (P1-P4) and pure cluster (Ml, 
Dl and Tl) geometries are used as reference systems 
throughout this work. A simple, representative structure 
for a polymer-cluster hybrid system has one Ag„ clus- 
ter attached to the polymer-monomer unit, for example, 
R-COO"Ag+ (M2, D2 and T2). We also consider more 
complicated structures where the Ag„ cluster is between 
the carboxylate groups of two such polymer-monomers 
(M7 and D6). Here, the Ag„ cluster is bound to four 
oxygen atoms, which create a stabilizing environment for 
it. 

The case of longer polymer chains bonded to Agi 
monomers and Ag 2 dimers is also addressed. Two ques- 
tions arise regarding the spectral properties of these sys- 
tems. First, does the length of the polymer chain itself 
have an effect on the absorption characteristics of the 
Ag„ cluster (M2 vs. M3 or M4, and D2 vs. D3 or D4)? 
Second, are the clusters attached to nearby branches of 
the system interacting with each other (M2-M4 vs. M5 
and M6 or M7 vs. M8, and D6 vs. D7)? 

When considering longer polymer chains, the concept 
of tacticity plays an important role in determining the 
lowest-energy structures. For a pure polymer, both iso- 
(P3) and syndiotactic (P4) isomers were considered. A 
similar comparison can be made between M3 and M4 
as well as between D3 and D4, even though the chain 
lengths are not equal. All of the structures with three 
polymer units in a row were chosen to be syndiotactic 
for two reasons. First, there is experimental evidence to 
suggest that the systems are likely to be syndiotactic 10 . 
Second, isotactic polymer chains with Ag„ clusters at- 
tached to neighboring polymer units tend to aggregate 
to form larger clusters. Thus, isotactic polymer chains 
do not favor the formation of stable hybrid structures 
with the small Ag„ clusters that are studied here. In the 
case of Ag2 clusters, the neighboring clusters aggregate 
even when they are attached to the next-neighbor poly- 
mer units for both iso- and syndiotactic polymers. The 
unstable structure D5 is shown in Fig. 1 as an example to 
high-light this point. In comparison, the structure D7 in 
which both dimers have stable anchorings to two polymer 
chains is found stable. 

The relative energies of the different structures were 
not studied in detail as it is not possible to compare 
the energetics of systems having different numbers and 
types of atoms. A direct energy comparison can only be 
made between the structures M2 and D6, and between 
the structures M5 and D7. In these cases, the D6 and 
D7 structures are favoured in energy by at least 1 eV per 
Ag atom. 



III. COMPUTATIONAL METHODS 

The SIESTA program u , which implements the 
standard Kohn-Sham self-consistent density-functional 
method, was used in spin-dependent calculations for the 
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D5 (Not stable) D6 Tl T2 



FIG. 1. (Color) Schematic pictures of studied structures (not to scale). Black, dark red, tiny white and small yellow dots 
represent PMAA units, their carboxylate groups, the hydrogen atoms of the carboxyl groups, and Ag atoms, respectively. 
The bond lengths or bond angles in the pictures do not represent the calculated values except for the detailed illustration of 
polymer-monomer PI, where large grey, large dark red and small light grey dots represent the carbon, oxygen and hydrogen 
atoms, respectively. 
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geometry optimizations. The Perdew-Burkc-Ernzcrhof 
(PBE) generalized gradient approximation (GGA) 42 , 
Troullier-Martins pseudopotcntials 1 , and the double 
zeta plus polarization (DZP) basis-set were also used in 
this work. The SIESTA program optimizes the geome- 
tries of the systems through a series of atomic displace- 
ments, until the calculated forces on the atoms are less 
than 0.01 eV/A. This is the standard procedure for geom- 
etry optimization using the SIESTA code, in which the 
numerical "noise", e.g., due to the real-space grid in the 
potential calculation, prevents the geometry optimiza- 
tion from being stuck at high-symmetry saddle points. 
Hence, in this work, Hessian analysis was not used for 
the structural optimization. 

The TDDFT program OCTOPUS 44 (version 
3.0.0/3.0.1) was applied in all further ground-state 
and excited-state calculations. Real-time TDDFT 
(RT-TDDFT) in its various flavors has been successfully 
used to determine the spectral properties of a large 
variety of both metal clusters and organic compounds 4 '. 
Calculating optical absorption spectra by TDDFT and 
comparing them with the measured spectra has proven 
to be an efficient tool for determining cluster and 
molecule structures, and even distinguishing between 
different isomers (see, e.g., Refs. 46 ' 47 ). In most cases, a 
good agreement with experiments up to a few tenths of 
an eV in spectral features is achieved. 

In this work, we used the RT-TDDFT method as 
implemented in OCTOPUS to calculate the absorption 
spectra in the linear response regime. The RT-TDDFT 
implementation follows Rcf . 1 , where the system is in- 
stantaneously excited by applying an electric held, and 
the occupied Kohn-Sham orbitals are then propagated 
in time. The photoabsorption cross-section follows from 
the induced dipole moments, and is calculated sepa- 
rately for each of the three spatial directions. The RT- 
TDDFT method should not be confused with linear- 
response TDDFT (LR-TDDFT) methods, such as the 
Casida method 11 ! . LR-TDDFT methods are superior for 
calculations of small systems. However, due to their poor 
scaling, the RT-TDDFT method becomes advantageous 
for large systems as it scales linearly with respect to the 
electron number. In this study, the Casida LR-TDDFT 
method within OCTOPUS was used to resolve the spe- 
cific transitions contributing to the absorption spectra 
for systems Ml, M2, Dl, and D2. 

The PBE exchange-correlation functional was used in 
the TDDFT calculations, except in the Casida calcu- 
lations, where the available kernel in OCTOPUS cor- 
responds to the local density approximation (LDA). 
Troullier-Martins pseudopotcntials 43 were used in all of 
the TDDFT calculations in this work. The calculations 
in the OCTOPUS program are based on a real-space 
grid. We used a spacing of 0.23 A for the real-space grid, 
and a simulation radius around each atom of 8 A, which 
are typical values in real-space RT-TDDFT calculations. 
The validity of the grid-spacing was verified by reference 
calculations for the test structure M2 using grid-spacings 



TABLE I. Properties of the studied systems. Even (odd) de- 
notes even (odd) number of electrons per Ag n cluster. All 
bond lengths are given in Angstroms. Cluster-cluster dis- 
tances are in italics. 
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2.33-2.34 




M6 


Even 


2.32-2.35 




M7 


Odd 


2.25-2.27 




M8 


Odd 


2.25-2.28 


6.62 


Dl 


Even 




2.581 


D2 


Odd 


2.27-2.32 


2.709 


D3 


Odd 


2.26-2.28 


2.710 


D4 


Odd 


2.27-2.28 


2.709 


D6 


Even 


2.17 


2.781 


D7 


Even 


2.17 


2.786, 6.52 


Tl 


Odd 




2x2.644, 3.023 


T2 


Even 


2.24-2.25 2x2.74, 2.723 



of 0.17 A and 0.20 A. All three grid -spacings produced 
similar absorption spectra up to an accuracy of about 
0.1 eV in the relevant low-energy region of about 0-7 eV. 
All relevant characteristics of the spectra were clearly re- 
solved using this accuracy, as the low-energy excitations 
of small metal clusters are well-spaced. The obtained 
numerical accuracy is also within the general accuracy of 
a few tenths of an eV for this method 23 ' 46 . The time- 
step for the calculation was 0.0025 h/eV and the number 
of time-steps was 20000 for the pure Ag„ clusters and 
the PMAA-monomer structure and 15000 for the rest of 
the structures, resulting in a total propagation time of 
50 ft/cV w 33 fs / 37.5 H/eV « 25 fs. In the Casida 
calculations, the number of unoccupied states that were 
taken into account was 4.2 times the number of occupied 
states for the largest structure D2 and at least 8 times 
the number of occupied states for the other structures. 
These numbers were considered adequate to produce reli- 
able results for the D2 structure and high-quality results 
for all of the other systems. 

IV. GROUND-STATE RESULTS 

The calculated bond lengths for the Ag containing 
structures are shown in Table I. In the case when the 
dimer is attached to the polymer, the Ag-Ag bonds are 
slightly lengthened (Dl vs. D2-D7). For the trimer, the 
Ag-Ag bond lengths in the triangle become almost equal 
(Tl vs. T2). The Ag-0 bond lengths vary only slightly 
between 2.17-2.42 A in the different structures, which is 
somewhat longer than the 2.04 A bond length in bulk 
silver oxide Ag20 . 

The ground-state calculations also give the HOMO- 
LUMO energy gaps, which can be considered as zcroth- 
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order approximations to the first excitation energy. We 
found that the structures of the same type, i.e., those 
having the same cluster size and the same number of 
polymer chains (for example, M2-M6 and D2-D4), have 
energy gaps of similar size. On the other hand, a 
strong gap variation exists between structures of differ- 
ent types. An odd number of electrons per Ag„ cluster is 
found to introduce a HOMO-LUMO energy gap smaller 
than 1.5 eV in all of the polymer-containing structures, 
whereas all polymer-containing structures with an even 
number of electrons per Ag„ cluster have gaps that are 
larger than 1.85 eV. Even though M8 has a total num- 
ber of electrons that is even, the odd number of elec- 
trons per Agi monomer in this system leads to local 
magnetic moments in the vicinity of the two Ag atoms, 
and to a HOMO-LUMO energy gap that is smaller than 
1.0 eV. For comparison, Refs. 51 ' 52 also report on cal- 
culated HOMO-LUMO energy gaps smaller than 1 eV 
for odd-electron Ag„ clusters within the size range of 
n = 2 — 20. The exception to this trend is the single 
silver atom (Ml in our notation), which has a gap that 
is greater than 3 eV. 



V. EXCITED-STATE RESULTS 

The absorption spectra for the pure polymer structures 
are given in Fig. 2. Each system shows a similar spectrum 
with an absorption edge above 4.5 eV. In the energy range 
of 4.5-7 eV, several peaks appear with slightly varying 
patterns depending on the type of structure. The spectra 
of structures P3 and P4, which differ only by tacticity, 
are almost identical in this energy range. The largest 
spectral strengths for these systems appear at higher en- 
ergies ranging up to several tens of eV's. This region, 
which is above the ionization energies of these systems, 
is not included in this figure as it is beyond the range of 
application of the TDDFT method. 

The absorption spectra corresponding to structures 
with monomers Agi, dimcrs Ag2 and trimcrs Ag3 
are shown in Figs. 3, 4 and 5, respectively. The 
pure Ag„ clusters, Ml, Dl and Tl, have their main 
lowest excitation energies at 3.8 eV, 3.0/4.5 eV and 
2.5/2.8/3.3/3.6 eV, respectively, in excellent agreement 
with previous studies 28 ' 29 . 

The inclusion of the PMAA monomer in structures M2 
and T2 causes only slight changes of a few tenths of an 
eV in the position of the lowest excitation and a slight de- 
crease in its intensity relative to those of the Ag monomer 
Ml and trimer Tl, respectively (cf. Ml and M2 in Fig. 3 
and Tl and T2 in Fig. 5). In addition, the spectral fea- 
tures of the polymer above 4.5 eV appear in both the M2 
and T2 cases. In contrast, for the dimer D2 in Fig. 4, the 
main low-energy excitations are perturbed with new ex- 
citations appearing at surprisingly low energies of 1.7 eV 
and 2.1 eV. From the individual directional components 
of the spectrum (not shown in the figure) , it can be con- 
cluded that these two low-energy excitations are caused 
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FIG. 2. (Color) Photoabsorption spectra of the pure polymer 
structures. The absorption strengths are scaled down by a 
factor of 1/2 for P2 and 1/3 for P3 and P4 to compensate for 
the different system sizes. 
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FIG. 3. (Color) Photoabsorption spectra of the structures 
with Agi monomers. The absorption strengths are scaled 
down by a factor of 1/2 or 1/3 for systems containing 2 or 3 
Ag atoms, respectively, to compensate for the different system 
sizes. Note the different y-axis scales. 
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FIG. 4. (Color) Photoabsorption spectra of the structures 
with Ag2 dimers. The absorption strength is scaled down by 
a factor of 1/2 for D7 to compensate for the different system 
size. 



by charge oscillations in the dimer-axis direction, whereas 
the excitations at 2.7-4.5 eV are along all three spatial 
directions. To compare, the pure dimer Dl excitation 
at 3.0 cV stems from the dimer-axis direction, and the 
doubly-intense excitation at 4.5 eV arises from the two 
other degenerate spatial directions perpendicular to the 
dimer-axis. Hence, we conclude that both low-energy 
peaks in the pure dimer spectrum have been destroyed 
and are replaced by a new pattern of multiple low-energy 
peaks in the spectrum of the D2 structure. We conjec- 
ture that this change in the low-energy spectrum arises 
due to the odd number of electrons in the vicinity of 
the Ag„ cluster corresponding to the perturbation in the 
electronic structure induced by the polymer. Such a per- 
turbation is already reflected in the ground-state results, 
for example, in the widths of the HOMO-LUMO energy 
gaps. 

Elongating the polymer chain causes only minor 
changes in the spectra. For example, compared to the 
single polymer-monomer systems M2 and D2, the spec- 
tra for the corresponding asymmetric structures M3 and 
D3 containing chains of two polymer-monomers, and the 
structures M4 and D4 containing chains of three polymer- 
monomers, differ only slightly (see Figs. 3 and 4). Neigh- 
boring Ag„ clusters do not significantly affect the spec- 
trum either. This can be seen, for example, in Fig. 3 
by comparing the spectrum of the M5 structure consist- 
ing of the syndiotactic chain of three polymer-monomers 



and the Ag atoms attached to next-neighbor polymer- 
monomers, with that of the structure M6 with the syn- 
diotactic chain of three polymer-monomers and the Ag 
atoms attached to each of the polymer-monomers. Only 
a slight redshift of about 0.1 eV and an extra shoulder 
at about 2.7-2.9 eV are observed in the latter case. To 
summarize, adding polymer units or even Ag atoms in 
M3-M6 does not change the fundamental trends in the 
spectra compared to M2. The same holds for the dimer 
structures D2-D4. 

When the Ag„ clusters are surrounded by two PMAA 
chains, dramatic changes in the absorption spectra are 
observed. For example, the low-energy excitations of M7 
and M8 correspond to the odd number of electrons per 
Ag n cluster in the hybrid structure. This is similar to the 
odd electron number effect seen in the D2, D3 and D4 
systems, which resulted in low-energy excitations. For 
the Ag monomer structure M7 (Fig. 3), the main low- 
energy excitation appears at 2.2 cV. Thus, this excita- 
tion is about 1.3 eV lower in energy than the low-energy 
excitations in M2-M6. The structure M8, which contains 
two such clusters at a mutual distance of 6.6 A, shows 
a moderate splitting and a slight redshift of the main 
excitation with respect to M7. 

In the case of dimers that are surrounded by two poly- 
mer chains (D6 and D7, Fig. 4), an opposite effect is 
observed. Namely, the absorption edge is brought up 
to over 4 eV. The low-energy excitation of the pure Ag 
dimer Dl (Fig. 4) at about 3 eV in the dimer-axis direc- 
tion is totally suppressed in the spectra for the D6 and 
D7 systems. The main excitations at about 5 eV and 
6.5 eV for the D6 and D7 systems arise from the direc- 
tion of the "bridges" between the two polymer chains. 
Thus, the presence of the polymer causes quenching of 
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the low-energy absorption of the Ag2 dimcr. Because no 
quenching was observed when Ag„ clusters were attached 
to single polymer chains (M2-M6, T2), we deduce that 
the formation of the symmetric structure of the Ag„ clus- 
ter and the four oxygen atoms in the D6 and D7 systems 
plays a key role in the suppression of certain directional 
vibration modes. 

To examine the origin of the shifts in the excita- 
tion spectra, the photoabsorption spectra of the small- 
est structures Ml, M2, Dl and D2 were calculated using 
the Casida method. This enables us to determine which 
orbitals are involved in the excitations, as such informa- 
tion is not directly available from the RT-TDDFT calcu- 
lations. The compositions of the main low-energy excita- 
tions of M2, Dl and D2 are shown in Figs. 6a-8a. For the 
pure Ag monomer Ml, the transition from the HOMO to 
the triply degenerate LUMO makes a 95% contribution 
to the low-energy excitation. Note that an almost con- 
stant shift of about 0.3 eV in transition energies between 
Figs. 3-4 and 6a-8a is caused by the use of the LDA ker- 
nel instead of the PBE kernel in the Casida calculations. 
To demonstrate the extent of the electronic perturbation 
induced by the polymer, the orbital illustrations of a few 
orbitals around the gap are also shown in Figs. 6b-8b, 
where the orbitals are illustrated as constant-value sur- 
faces of the electron wavefunctions. 

Fig. 6a shows that the main low-energy excitation in 
the absorption spectrum of M2 corresponds to the tran- 
sitions from the HOMO-1 and HOMO-2 orbitals to the 
LUMO orbital. As shown in Fig. 6b, these orbitals are 
localized over the carboxylate group and the Ag atom, 
and even over the carbon skeleton. This differs clearly 
from the pure Ag atom HOMO-LUMO excitation from 
the s-orbital to the triply degenerate p-orbital. 

For the Ag dimer Dl, the LUMO is an anti-bonding 
orbital and the degenerate LUMO+1 and LUMO+2 arc 
7r-bonding orbitals as can be seen in Fig 7. A similar 
pattern with a clear interaction between the Ag dimcr 
and the carboxylate group is seen in the orbitals of the 
D2 structure in Fig. 8. The difference between Dl and 
D2 is that D2 has an odd number of electrons, and a net 
spin resulting in an absorption spectrum with several ex- 
citations at lower energies. Previous results by Mitric et 
al. 35 for Ag n cluster-tryptophan structures show similar 
features to our findings. For the Ag dimer attached to 
the carboxylate group of a tryptophan molecule, an exci- 
tation at an anomalously low energy of about 2.3 eV was 
observed'' 1 . This effect was associated with a HOMO- 
LUMO transition that is similar to the one described 
above for our system, D2. 

Next, we consider the quenching of the lowest exci- 
tation peak when moving from the dimer Dl to struc- 
tures D6 and D7 (See Fig. 4). The increase in the bond 
length between the Ag atoms is far too small to explain 
the quenching. Instead, we note that the HOMO-LUMO 
transition makes a 94% contribution (Fig. 7) to this ex- 
citation. The quenching that is observed is therefore due 
to the destruction of the a shape of the HOMO, which 



happens already in the D2 structure as seen in Fig. 8. 

Finally, we want to consider the charge transfer as a 
possible source of error in our calculations. The effect of 
charge transfer can lead to errors of several eV's in exci- 
tation energies for the commonly used (semi) local func- 
tionals such as PBE oi . The hybrid functionals such as 
B3LYP do not perform much better. To probe the relia- 
bility of the PBE functional in our calculations, we have 
used the criterion by Peach et al. to estimate the amount 
of charge transfer . In this criterion, Peach et al. define 
the spatial overlap A as attaining values between and 
1, so that large values correspond to low charge transfer. 
According to Peach et al. the PBE results are predicted 
to be reliable for A > 0.4. We found that for the struc- 
ture M2, A is for the lowest energy excitation within the 
range of 0.50-0.54. For the structure D2, the A values 
for the five low-energy excitations shown in Fig. 8a are, 
in the order of increasing excitation energy, within the 
ranges of 0.71-0.73, 0.60-0.63, 0.60-0.62, 0.43-0.46 and 
0.61-0.64, respectively. The uncertainty in these values 
is due to the inclusion of a limited number of orbitals in 
the calculation of A. From these values and the orbital 
illustrations in Figs. 6 and 8, we conclude that the char- 
acter of the transitions is mainly local, and that the PBE 
functional can be safely used to study the main features 
of the excitations in these systems. To see whether the 
same conclusion can be drawn for our larger systems, we 
have examined the electronic structure of M5. We ob- 
serve that the orbitals near the gap of this system have 
shapes similar to those of M2 (See Fig. 6). Moreover, 
each of the orbitals of M5 is totally or mainly localized 
in the vicinity of one Ag atom only. Supported by the 
previous result that the photoabsorption spectra changed 
only slightly with the chain length, we can conclude that 
the excitations of M5 are also local in character. Using 
the structure D2 as a reference, a similar analysis was 
also performed for D6, which yielded the same result. As 
the rest of the structures have a strong resemblance to 
these test cases, we do not expect the charge transfer to 
be a major source of error in any of the excitations in the 
studied structures. 



VI. CONCLUSIONS 

We have studied how the photoabsorption spectra of 
Agi , Ag2 and Ag3 clusters are affected by the presence of 
the PMAA polymer. The structure optimization shows 
that there exists stable configurations in which the Ag„ 
clusters are bound to the carboxylate group of the poly- 
mer. 

We observe that the bonding of the polymer to the Ag„ 
clusters has an extensive effect on the absorption char- 
acteristics in certain cases. Namely, when the number of 
electrons per Ag„ cluster is odd, the electronic structure 
and thereby the spectral features of the Ag„ clusters are 
strongly perturbed. For example, for the PMAA _ Ag2~ 
systems with an odd number of electrons, the minimum 
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FIG. 6. (Color) (a) Energy level diagram of the main low- 
energy excitation at energy E and (b) orbital illustrations of 
the structure M2. Light blue represents positive and dark 
blue negative parts of the wavef unction. 



absorption energy is lowered by almost 1.5 eV resulting 
in absorption in the visible range. Suppression of the 
lowest excitations is observed when an Ag ra cluster with 
an even number of electrons, such as Ag 2 , is surrounded 
by two carboxyl groups, so that the cluster is bonded to 
four oxygen atoms. 

Our results reveal that the absorption features of the 
smallest Ag„ clusters are very sensitive to their chemical 
environment. This suggests that these clusters are viable 
candidates for materials with tailored optical properties. 
Already, there exists experimental evidence on the ability 
to tune the absorption and emission spectra of small Ag„ 
cluster-polymer hybrid systems. The results presented 
here give theoretical confirmation of those findings, thus 
showing such systems to be promising candidates for de- 
veloping functional nanomaterials with adjustable prop- 
erties. 



© C LUMO+(1-2) 




FIG. 7. (Color) (a) Energy level diagram of the main low- 
energy excitations at energies E and (b) orbital illustrations 
of Ag dimer Dl. The excitation at E = 4.83 eV is doubly 
degenerate as the LUMO+1 and LUMO+2 states are degen- 
erate orbitals. Light blue represents positive and dark blue 
negative parts of the wavefunction. Reference energy level 
diagrams can be found in Refs. ' . 
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